An efficient edge based data structure for the compressible Reynolds-averaged Navier–Stokes equations on hybrid unstructured meshes

An efficient edge based data structure has been developed in order to implement an unstructured vertex based finite volume algorithm for the Reynolds-averaged Navier–Stokes equations on hybrid meshes. In the present approach, the data structure is tailored to meet the requirements of the vertex based algorithm by considering data access patterns and cache efficiency. The required data are packed and allocated in a way that they are close to each other in the physical memory. Therefore, the proposed data structure increases cache performance and improves computation time. As a result, the explicit flow solver indicates a significant speed up compared to other open-source solvers in terms of CPU time. A fully implicit version has also been implemented based on the PETSc library in order to improve the robustness of the algorithm. The resulting algebraic equations due to the compressible Navier–Stokes and the one equation Spalart–Allmaras turbulence equations are solved in a monolithic manner using the restricted additive Schwarz preconditioner combined with the FGMRES Krylov subspace algorithm. In order to further improve the computational accuracy, the multiscale metric based anisotropic mesh refinement library PyAMG is used for mesh adaptation. The numerical algorithm is validated for the classical benchmark problems such as the transonic turbulent flow around a supercritical RAE2822 airfoil and DLR-F6 wing-body-nacelle-pylon configuration. The efficiency of the data structure is demonstrated by achieving up to an order of magnitude speed up in CPU times.     

A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

A survey of surface tension,molar volume and density for Sn–Ag–Cu–Bi–Sb quinary alloys as lead-free solders

ABSTRACT

By applying the geometric models and the theoretical equation, the surface tension, the molar volume and the density were studied. The empirical calculations were carried out in temperature range 623?K?≤?T?≤?1123?K. Only few thermophysical properties were estimated for eight quinary alloys: Sn3.55Ag0.5Cu3Bi3Sb, Sn3.48Ag0.5Cu3Bi5Sb, Sn3.48 Ag0.5Cu5Bi3Sb, Sn3.40 Ag0.5Cu5Bi5Sb, Sn3.53Ag1Cu3Bi3Sb, Sn3.46Ag1Cu3Bi5Sb, Sn3.46Ag1Cu5Bi3Sb, Sn3.38Ag1Cu5Bi5Sb. The results show that surface tension and density have a linear appearance for all temperatures. We have also studied the influence of the composition and temperature in the studied alloys. The obtained theoretical results are compared with the experimental ones and with the conventional Pb–Sn welds.     

Physical aging in amorphous poly(ethylene furanoate): Enthalpic recovery,density, and oxygen transport considerations

The current work utilizes three separate techniques to study the physical aging process in amorphous poly(ethylene furanoate) (PEF), which is a recently introduced engineering thermoplastic with enhanced properties compared to petroleum‐sourced poly(ethylene terephthalate). Differential scanning calorimetry aging experiments were conducted at multiple aging temperatures and times, and the resultant enthalpic recovery values compared to the theoretical maximum enthalpy loss evaluated from calculations involving extrapolation of the equilibrium liquid line. Density measurements reveal densification of the matrix for the aged versus unaged samples, and provide an estimate for the reduction in free volume for the aged samples. Complementary oxygen permeation and pressure‐decay sorption experiments provide independent verification of the free volume reduction mechanism for physical aging in glassy polymers. The current work provides the first detailed aging study for PEF. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53, 389–399     

T1Effects in Sequential Dynamic Susceptibility Contrast Experiments

Residual effects of an initial bolus of gadolinium contrast agent have been previously demonstrated in sequential dynamic susceptibility contrast MR experiments. While these residual effects quickly reach a saturation steady state, their etiology is uncertain, and they can lead to spurious estimates of hemodynamic parameters in activation experiments. The possible influence ofT₁effects is now investigated with experiments in whichT₁weighting is varied as well as with serial regionalT₁measurements. Little evidence for significant residualT₁effects is found, suggesting instead that susceptibility effects underlie these observations. An initial saturation dose of contrast agent minimizes this effect.     

Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

乌桕类可可脂结晶过程中表现体积变化

为了解乌桕类可可脂(CTCBE)缓冷结晶横向胀罐爆裂原因，研究CTCBE结晶过程中的膨胀特性，采用流体静力法测试了不同等温结晶条件下CTCBE结晶形成的表观密度和表观体积及其横纵向膨胀状况。结果表明，CTCBE缓冷结晶后形成外部和中部两个晶区。在5-25℃的各等温结晶范围内，当结晶温度增加时，外部结晶区域减小、表观密度变化不大；中部结晶区域增大、表观密度明显减小，CTCBE的表观体积、横向、纵向膨胀率均增大，膨胀横向大于纵向。进一步表明CTCBE在自然缓冷固化的表观体积根本上由其结晶温度决定。为有效防止横向胀罐，自然缓冷固化温度应低于10℃。     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

Statistical thermodynamics evaluation of polymer–polymer miscibility

The Simha and Somcynsky (S–S) statistical thermodynamics theory was used to compute the solubility parameters as a function of temperature and pressure [δ = δ(T, P)], for a series of polymer melts. The characteristic scaling parameters required for this task, P*, T*, and V*, were extracted from the pressure–temperature–volume (PVT) data. To determine the potential polymer–polymer miscibility, the dependence of δ versus T (at ambient pressure) was computed for 17 polymers. Close proximity of the δ versus T curves for four miscible polymer pairs: PPE/PS, PS/PVME, and PC/PMMA signaled the usefulness of this approach. It is noteworthy, that the tabulated solubility parameters (derived from the solution data under ambient conditions) propounded the immiscibility of the PVC/PVAc pair. The computed values of δ also suggested miscibility for polymer pairs of unknown miscibility, namely PPE/PVC, PPE/PVAc, and PET/PSF. In recognizing the limitations of the solubility parameter approach (the omission of several thermodynamic contributions), these preliminary results are auspicious because they indicate a new route for estimating the miscibility of any polymeric material at a given temperature and pressure. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 2909–2915, 2004